Measurements of Deuteron Photodisintegration up to 4.0 GeV 
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The first measurements of the differential cross section for the d{^,p)n reaction up to 4.0 GeV 
were performed at Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. We 
report the cross sections at the proton center-of-mass angles of 36°, 52°, 69° and 89°. These results 
are in reasonable agreement with previous measurements at lower energy. The 89° and 69° data show 
constituent-counting-rule behavior up to 4.0 GeV photon energy. The 36° and 52° data disagree 
with the counting rule behavior. The quantum chromodynamics (QCD) model of nuclear reactions 
involving reduced amplitudes disagrees with the present data. 
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To reconcile low energy and high energy descriptions of 
hadronic matter, nuclear physics must determine when it 
is justified to make a transition from meson-nucleon de- 
grees of freedom to quark-gluon degrees of freedom in the 
description of a nuclear reaction. The QCD content of 
nuclei was studied first by Brodsky and Chertok 0. A 
possible signature for this transition is that the reaction 
cross section begins to scalej^ at some incident energy. If 
scaling were indeed observed, characterization of the ap- 
proach to scaling would be essential to understand how 
the dynamics are simplified. High energy two-body pho- 
todisintegration of the deuteron (-fd — > pn) is particu- 
larly well suited for these studies because it is amenable 
to theoretical calculation and relatively high momentum 
transfer to the constituents can be achieved at relatively 
modest photon energies Q . 

Previous measurements |^ for the d( , y,p)n reaction 
indicate the onset of scaling behavior at a photon energy 
of 1 GcV at a reaction angle of 9 cm = 90°; however, this 
limited data set does not show scaling at other reaction 
angles. Measurement of the angular distributions of scal- 
ing thresholds is an important part of characterizing the 
reaction process. In this Letter, we present new data to 
address this issue. 

Predictions for the energy dependence of the cross sec- 
tion at high energies are given by constituent counting 
rules H , the reduced nuclear amplitude analysis (RNA) 
P|, the asymptotic meson exchange model (AMEC) |7]] 
and the quark-gluon string model (QGS) Q. Constituent 
counting rules predict that the energy dependence for the 
two-body exclusive reaction cross section should be given 
by 



da/dt = 



h(9 cm ) 



(1) 



where the Mandelstam variables s and t are the square of 
the total energy in the center-of-mass frame and the mo- 
mentum transfer squared in the s channel respectively. 
The symbol n denotes the total number of elementary 
fields in the initial and final states, and n — 2 is 11 for the 
d(7,p)n reaction. The quantity h(9 cm ) depends on de- 
tails of the dynamics of the process. Constituent count- 
ing rules are believed valid at energies much greater than 
the masses of the participating particles, perhaps in the 
perturbative QCD region. The previous data at 90° scale 
according to these rules with a photon energy threshold 
of only 1 GeV. 

In the RNA approach, the amplitude is described in 
terms of parton exchange between the two nucleons. The 
low energy components responsible for quark binding 



1 Scaling in this context implies a dependence on a reduced 
set of kinematic variables indicating a simplification of the 
reaction dynamics. 



within the nucleons are removed by dividing out the 
empirical nucleon form factors. While the reduced nu- 
clear amplitude analysis appears to describe the electron- 
deuteron elastic scattering cross section above a momen- 
tum transfer of 1 GeV/c B, it does not give a good 
description of the previous deuteron photodisintegration 
data. This is surprising because this model is expected to 
approach scaling at lower energies than the constituent 
counting rules. The QGS model is based on Regge phe- 
nomenology and is expected to be valid at small values 
of t where the parameters are best determined. Thus, 
the QGS model is expected to be valid at small reaction 
angles, while most of the existing data is at large angles. 

The traditional meson exchange models [p|Jl0[| describe 
the data at energies below 1 GeV, but are problematic 
above 1 GeV. However, the AMEC departs from the 
conventional approach in that an asymptotic description 
of the nucleon- nucleon interaction is used [[llj . Although 
this model appears to predict the observed energy depen- 
dence at 9 cm = 90°, it cannot yet reproduce the magni- 
tude of the cross section. 

In this Letter we present new results from the d{pi,p)n 
reaction at proton center-of-mass angles of 36°, 52°, 69° 
and 89°. These data overlap with existing SLAC mea- 
surements, extending them to higher energy and provid- 
ing more complete angular coverage. 

Because the cross sections are much less than 1 nb/sr, 
it was essential to have a high-current, high duty fac- 
tor electron beam of multi-GeV energy, and to use well- 
shielded spectrometers of large solid angle. This was 
achieved by performing the experiment in Hall C at 
Thomas Jefferson National Accelerator Facility (Jeffer- 
son Lab). A 25 uA continuous (CW) electron beam in the 
energy range 0.845 to 4.045 GeV in steps of 0.8 GeV was 
incident on a 6% copper radiator to create an untagged 
photon beam. The resulting electron and bremsstrahlung 
beam impinged on a 15-cm liquid deuterium target. Be- 
cause the 7G? — > pn reaction is a two-body process, the 
photon energy could be reconstructed from the mea- 
sured final-state proton momentum and scattering angle. 
Events with pion production were excluded by accepting 
only the protons with the highest momenta. Protons 
were detected with the High Momentum Spectrometer 
(HMS) @. 

The HMS had a solid angle of 6.8 msr and a momen- 
tum acceptance of 18%. The detector system consisted 
of a plastic scintillator hodoscope, drift chambers, and a 
gas Cerenkov counter. The hodoscope was used to form a 
trigger and also to provide the time-of-flight information 
for particle identification. Drift chambers were employed 
to measure particle trajectories in order to calculate the 
momenta and reaction angles. It is important to distin- 
guish protons from pions and deuterons. A gas Cerenkov 
counter, filled with C4-F10 gas, was used in combination 
with time-of-flight (TOF) information to separate pions 
and protons above 2.8 GeV/c, while TOF alone was used 
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for pion rejection at lower momenta. Deuteron rejection 
was accomplished by a TOF cut at all momenta. 

Background contributions from the target windows 
were removed by placing cuts on the reconstructed target 
position and subtracting the yield obtained with a cell 
of identical dimensions that was either empty or filled 
with liquid hydrogen to simulate bremsstrahlung in the 
deuterium. This target length cut led to a 3% error for 
the target thickness. Deuterium and empty target data 
were taken alternately during the experiment. Data were 
taken with the hydrogen target at some of the kinematic 
settings to cross check the procedure of the empty tar- 
get subtraction. The yield from electrodisintegration was 
measured by repeating the procedure without the radi- 
ator present. This yield was treated as a background 
and was subtracted from the photodisintcgration yield 
with an energy-dependent correction factor to take into 
account the modification of the electron beam flux and 
energy distribution by the radiator fl3| . An additional 
systematic error of 3% was assigned to this background 
subtraction procedure. 

The photon energy bin limits were chosen to kinemat- 
ically eliminate protons from photopion production pro- 
cesses, and to eliminate the bremsstrahlung end point, for 
which the photon flux is less well known. This photon en- 
ergy bin varied in width from 60 MeV at 36° to 115 MeV 
at 89°. The bremsstrahlung photon flux was calculated 
with an estimated 3% uncertainty using the thick-target 
bremsstrahlung computer code developed by Belz Q , 
which was cross-checked by using the codes of Matthews 
and Owens . The spectrometer solid angle for the ex- 
tended target was studied by comparing a Monte Carlo 
simulation with measurements made with both a mov- 
able solid carbon target and the extended cryotargets. 
In addition, the dependence of the results on the choice 
of optics parameters for the HMS was studied. An over- 
all error of 7% in the HMS acceptance was determined 
from these studies. A proton absorption correction was 
applied to compensate for the scattering in the spectrom- 
eter windows and the detector stack. This proton attenu- 
ation was measured [Jl5| to be (5.5 ±2.0) % by comparing 
singles H(e,e') and coincidence H(e,e'p) measurements. 
Corrections were also applied for the computer deadtime 
and the tracking efficiency. 

The overall systematic uncertainty is found to be < 
11.5%. The uncertainties from the beam current mea- 
surement, beam energy determination and photon en- 
ergy reconstruction in the measured quantity s xl dajdt 
were less than 3%. The uncertainty from the particle 
identification is < 5%. An additional systematic error of 
3% resulted from using a reduced solid angle in the HMS 
to avoid an obstruction from an HMS vacuum valve that 
was inadvertantly partially closed |16| , |l7j| . Thus, only 
half of the HMS solid angle was used in the analysis of 
the data. A separate experiment without the obstruc- 
tion verified that there was good agreement between the 



half acceptance of the experiment with the obstruction 
and the corresponding half acceptance for the experiment 
without the obstruction. A Monte Carlo simulation was 
used to determine the correction for "in-scattered" pro- 
tons from the valve into the open half of the spectrometer 
solid angle. The amount of this correction is less than 2% 
in the worst case. This selection of half the spectrometer 
solid angle resulted in a reduction of the observed reac- 
tion center of mass angle by approximately 1° from the 
central spectrometer setting. 

Figure 1 shows the world data of s 11 da/dt above ~ 0.4 
GeV for the jd — > pn reaction at 8 cm — 36°, 52°, 69° and 
89° as a function of the photon energy. The 89° data 
are shown in the top panel of Fig. 1. The SLAC NE17 
data Q exhibit scaling behavior starting at photon en- 
ergies around 1 GeV, corresponding to proton transverse 
momenta (pr) of ~ 1.0 GeV/c. For the d(fy,p)n reac- 
tion, pt 2 can be expressed as \E 1 MdSin 2 {9 cm ) with Md 
being the deuteron mass. The Jefferson Lab data are in 
good agreement with the previous SLAC data. The new 
data continue to show the scaling behavior up to the 
highest photon energy of 4.0 GeV. The differential cross 
section da/dt from the present work at photon energies 
above 1.0 GeV were fitted with the form of A/s n ~ 2 , and 
n — 2 = 11.1±0.3 was obtained from the fit. It is sur- 
prising that the counting rule appears to work so well 
considering the fact that the momentum transfer to the 
individual quark is below 1.0 (GeV/c) 2 for these mea- 
surements, where the strong coupling constant still varies 
significantly as a function of the momentum transfer. 

The solid line in Fig. 1 represents the meson-exchange 
calculation of Lee S, which is a standard calculation 
that reproduces the measured NN phase shifts up to 2.0 
GeV and is also constrained by photomeson jl9| produc- 
tion data. The calculation gives a reasonable description 
of the data below 0.5 GeV, but deviates above 1.0 GeV. 
The meson-exchange model calculation of Laget JToJ (not 
shown) reproduces the experimental data fairly well from 
the threshold up to about 0.6 GeV. The long-dashed line 
represents a reduced nuclear amplitude analysis || with 
a normalization factor chosen to agree with the data at 
E y = 1.6 GeV. This RNA curve falls below the high 
E 1 data and does not reach an asymptotic limit at these 
energies. The dotted line is a calculation (AMEC) by 
Nagornyi et al. H] with the normalization performed at 
1.0 GeV. Although this calculation gives a different en- 
ergy dependence to the cross section than that of the 
constituent counting rule, it does give an energy depen- 
dence very close to 1/s 11 in the energy region of 1.5 to 
3.5 GeV. More calculations of this type are necessary at 
the other angles. The QGS calculation j^] is shown as 
the dash-dotted line. It deviates from the data at large 
angles as expected. 

The data at 52° and 69° are shown in the center pan- 
els of Fig. 1. The 69° data appear to scale according to 
the constituent counting rules. A fit to this cross section 
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gives n — 2 — 10.9 ± 0.2 in good agreement with that at 
89°. At 52° the present data are in reasonable agreement 
with the previous SLAC measurement and have signifi- 
cantly improved statistical errors. The results deviate 
from the scaling behavior predicted by counting rules. A 
fit of the scaling formula to the cross section data gives 
n - 2 = 9.5±0.2. The QGS calculation (§ is in rea- 
sonable agreement with the data in the photon energy 
region between 1.0 to 4.0 GeV. 

The bottom panel in Fig. 1 shows the data at 9 cm = 
36°. Fitting the scaling formula to the Jefferson Lab 
cross section data above a photon energy of 1 GeV gives 
n - 2 = 9.6 ± 0.1. The highest p T from NE17 mea- 
surements at this 8 cm is 0.7 GeV/c. The experimental 
uncertainties and the limited energy region covered by 
the NE17 experiment preclude any conclusion with re- 
gard to the scaling behavior in the photon energy range 
below 3 GeV. The Jefferson Lab data are in reasonable 
agreement with the SLAC data |Q] when the large angu- 
lar dependence at forward angles is taken into account. 
The SLAC data were taken at 36.7°, while the present 
data were recorded from 35.4° at the highest energy to 
36.2° at the lowest energy. The previous data shown in 
Fig. 1 are not corrected for this difference in angles. 

Although a photon energy of 3.1 GeV at 36° corre- 
sponds to the same pt where the onset of s -11 scaling is 
observed at 90°, there is no evidence from the 36° data 
that scaling has set in (see arrows in Fig. 1). The RNA ||] 
deviates significantly from the data at 36°, but the QGS 
calculation Q gives a reasonable energy dependence. 

To summarize, the data at 9 cm = 89 — 90° continue 
to show scaling behavior up to 4 GeV. The first evi- 
dence is observed for a similar scaling behavior at 69°. 
The asymptotic meson exchange model is in fair agree- 
ment with the 90° data in the photon energy region of 
1.0 to 3.0 GeV. Thus far, no constituent counting scal- 
ing behavior is observed for the 36° or 52° data. The 
new data at 36° and 52° rule out a scaling threshold 
of Pt — 1-0 (GeV/c) 2 . The reduced nuclear amplitude 
analysis does not agree with the present data. The QGS 
model has an energy dependence that is in reasonable 
agreement with the present data at 36° and 52°. Further 
measurements at higher energies at the forward angles 
will test whether scaling has really set in and whether px 
is the physical quantity governing the onset of the scal- 
ing behavior. Experimental determination of the physical 
observable associated with the onset of scaling is essen- 
tial for understanding the underlying mechanism for the 
scaling behavior in the exclusive two-body process. 
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FIG. 1. s da/dt vs. _B 7 . The present work is shown as 
solid circles with statistical uncertainties only, the SLAC 
NE17 data are shown as open circles, the SLAC NE8 data 
are shown as open triangles, and other existing data Jl8| are 
shown as crosses. The solid line is the meson- exchange model 
calculation of Lee Jsjj . The long-dashed line is the RNA anal- 
ysis Q, and the dotted line is Nagornyi's j?j asymptotic me- 
son-exchange calculation. The dash-dotted line is the QGS 
calculation |8| . The arrows indicate the photon energies where 
Pt 2 = 1.0 (GeV/c) 2 . The previous data are shown above at 
nominal center of mass angles of 37° , 53° , and 90° . 
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